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1. Introduction 

In a recent paper we reported on the calculation and analysis of low-temperature 
series for the square lattice spin-1 Ising model. The series was derived to 79th order 
(in the variable u = exp[— J/fc^T]) using the finite-lattice method and employed 
a new algorithm which removed much of the memory-size restrictions of previous 
implementations. In this paper we report on a further extension of these series to 
order 113 for the specific heat and spontaneous magnetisation and order 112 for the 
susceptibility. The extension is obtained by direct calculation of the series to order 
99 and use of a new extrapolation procedure to extend the series by an additional 13 
or 14 terms. The improvement in the direct series derivation is due to a more efficient 
implementation of the algorithm and the use of parallel computation. The extrapolation 
procedure is similar to and inspired by work on directed percolation 0. 



2. The series expansion technique 



The Hamiltonian defining the spin-1 Ising model in a magnetic field h can be written: 
n = Jj2i^-a.a,) + hJ2{l-a.) (1) 

where the spin variable CTj = 0, ±1. The first sum is over nearest neighbour pairs and 
the second sum is over sites. The constants are chosen so the ground state (cXi = +1 V i) 
has zero energy. The low temperature expansion is based on perturbations from the 
fully aligned ground state. The expansion is expressed in terms of the temperature 
variable u = exp(— /3J) and the field variable /i = exp(— where P = l/ksT. The 
expansion of the partition function in powers of u may be expressed as 



Z=J2 



(2) 



n=0 



where '^nifJ') are polynomials in fi. It is more convenient to express the field dependence 
in terms of the variable x = 1 — fi and truncate the expansion at x"^ 

Z = Zo{u)+xZ,{u)+x^Z2{u) + (3) 

where Zn{u) is a series in u formed by collecting all terms in the expansion of Z 
containing factors of a;". Standard definitions yield the spontaneous magnetisation 



M{u) = M(0) + 



1 d\nZ 



(3 dh 



Zi{u)/Zo{u), 



h=0 



since x = in zero field. For the zero-field susceptibility we find 



dM 



dh 



d 



1 dZ 



h=0 



dh \(3Z dh 



/3 



h=0 



, Z2{u) _ Zi{u) _ / Zi{uY 
' Zo{u) Zo{u) \Zo{u)^ 



(4) 



(5) 
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The specific heat series is derived from the zero field partition function (via the 
internal energy U = —{d/d(5) InZg), 

dJJ rP' / d \ ^ 

CM = gf = /^'^l^^o = if^J)' U^j ln^o(«). (6) 

So in order to obtain the series expansion of the specific heat, spontaneous 
magnetisation and susceptibility it suffices to calculate the three quantities Zo, Z\ and 

On the square lattice the infinite lattice partition function Z can be approximated 
by a product of partition functions Z^n on finite [m x n) lattices, 

Z{u) ~ JJ^ Zmn(w)"'"" with m <n and m + n < r, (7) 

m,n 

where r is a cut-off which limits the size of the rectangles considered. The weights a^n 
are small integers and are known explicitly for the square lattice. 



1 


if m + n = 


r 




-3 


if m + n = 


r 


- 1 


3 


if m + n = 


r 


- 2 


-1 


if m + = 


r 


-3 





otherwise 







Due to the symmetry of the square lattice one obviously has that Zmn = Znm so one need 
only consider the case m < n and change the weights amn appropriately, i.e., multiply 
by 2 if m < n. 

For the low-temperature expansion of the Ising model Zmn is calculated as the sum 
over all spin configurations on the finite lattice. All spins outside the mxn range are 
fixed at +1. The number of terms derived correctly with the finite lattice method is 
given by the power of the lowest-order connected graph not contained in any of the 
rectangles considered, which in this case are chains of sites all in the '0' state. From the 
Ising Hamiltonian we see that such chains give rise to terms of order 3r + 1 . For a given 
value of r the series expansion is thus correct to order 3r. 

The efficient way of calculating Zmn is by transfer matrix techniques. We refer to 
|l|] for a detailed description of the algorithm. For this work we used a more efficient 
implementation of the algorithm and a parallel computer and were able to derive the 
series directly up to a maximal cut-off r.m = 33. 

3. Extrapolation of series 

The series can be extended significantly via an extrapolation method similar to that of 
0. Consider the series for Zn{u). For each r < we use the finite- lattice method 



to calculate the polynomials Z^^riu) = J2j=o ^n,j,rU^ correct to As already 

noted these polynomials agree with the series for Z„ to 0{u^^'). Next, we look at the 
sequences dn,r,s obtained from the difference between successive polynomials 



U 



3r+l 



s _ ^3r+l 



rU . 



(9) 



s>0 



s>0 



The first of these correction terms dnfl,r is often a simple sequence which one can readily 
identify. In the case of Zq we find the sequence 



0,0,r 



1,2,6,18,52,138,338,778,1712, 



from which we conjecture 

do,o,r = -2*^+^ + {r^ + 3r^ + 2r + 18)/3, r > 2. 



(10) 



The formula for (io,o,r holds for all the r^ — 1 values that we calculated and we are very 
confident that it is correct for all values of r. As was the case in the higher-order 
correction terms dn,s,r can be expressed as rational functions of dnfl,r- Due to the form 
of the first correction term this leads to the general extrapolation formulae 

-| s+n+3 -j^ s+n 



dr. 



E 



QsW. 



+ ^ E Ks,jr'2\ r>s + 2. 
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slnl 



j=0 



The factors in front of the sums have been chosen so as to make the leading coefficients 
particularly simple. We were able to find formulae for all correction terms up to 5 = 13 
for Zq and Zi and up to 5 = 12 for The coefficients in the extrapolation formulae 
are listed in tables [l|-^. 

Table 1. Coefficients ao,sj and &o,sj in the extrapolation formula (O) for Zq. 
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22872 
-15300 
3708 



e337929e 4055122080 

-33706348 - 1962612264 

6885452 372299004 

-692008 -35838744 

39698 2015962 



288647444160 22547966414400 1918636884702720 

-128643760656 -9349955599536 -746245210622400 

22859538096 1569303632880 119023740391728 
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163920 
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- 18208734 
327676 
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-13364 
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5670 
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-8136639 
1736096j 
-46873 
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-23158S134427|- 
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-35390845} 
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-321585 
-14080 
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-734477218900425 
-406644S6659736i 
-6538249560737i 
1573709594180^ 
-73717696634^ 
13227448244li 
-8870615^ 
-3532691^ 
-300061751 



-31757793382715400 
-1919833204276777^ 
- 170426378459944 
32784296682271} 
687212762691^1 



13285766431>5| 
-4219620181.^ 
213929108^ 
34494495j 
-1097239 
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Table 2. Coefficients ai^s,j and bi^s,j in the extrapolation formula (|ll| ) for Zi. 



s/j 1 2 3 




5 


6 


7 


8 




10 


11 






-24 -1008 -76272 -5684112 


-457173792 


- 39103403040 


-3574330130880 


-348267251263680 




-3973220722801612800 


-462907224587766067200 


-56960092350846168422400 


-73848939500620523 79648000 


1006557478598701079008368400 


1 -10 242 26776 2132988 




14031817536 


1234927394784 


115716608013024 


11544174716755968 


1224242942187936384 


137714989912773327360 


16393639864515714443520 


2060204208928428438812160 


272692487745604926720576000 


2 -13 161 3890 -11104 


-10362700 


-1228736780 


-124161954720 


-12398829627432 


-1273546890227328 


- 136597570576055328 


-15386179110821861184 


-18233261SS928804380864 


-227313813910382067927168 


-29786106025581577422698880 


3 -5 -76 -2693 -88837 


-3049S40 


-100973728 


-1890947308 


191723836836 


40750032405880 


5702911364579304 


735898260522888000 


94356730717840349328 


12349621039021303100608 


1669025207886453938464336 


4-2 7 31? 15958 


526035 


23570715 


1144281164 


57449415410 


2729697798164 


90106090294292 


- 4306747906782568 


-1641298736512799072 


-300489154903171705008 


-48094001062232864660800 


5 2 -11 -793 


-35351 


-1674878 


-88419027 


-5114559137 


-315077153986 


-20142410271322 


- 1282219400735740 


-73711506188062096 


-2454291894283312912 


302424204333008932640 




IBll 


71730 


3837282 


232243641 


15355331735 


1081711840093 






475590666188289368 




7 2 


-25 


-2852 


-121842 


-7085208 


-490136073 


-36372405428 


-2855093133545 


-236164037257269 


-20420710128454372 


-1819694197681939020 


S 




35 


4536 




11069628 


891555753 


73516402976 


6324527102839 


572875721359139 


54403254827920427 


9 






-47 


-6917 


-240352 


- 14817384 


-1460017005 


-132977895083 


-12335475680011 


- 1208966177416998 


10 






-2 


61 


gss3 


286375 


16521220 


2223262547 


222732548753 


21808039004180 


11 








! 


-77 


-13556 


-300155 


-13327961 


-3238529375 


-355244044168 














95 


17932 




1130371 


4637002786 


13 












2 


-115 


-23241 


- 127937 


25630356 



11367045 
659431i 
130525i 

-32166^ 
56471 



374646825 
6912974; 
2055016; 
-703727^ 
-6373^ 
-4700^ 



13414424940 
- 6387524: 
70181556^ 



521931313260 
-9516059361 
2443225632} 
-104301900 
-49051827i 
3508470 



-4278680737j 
-1421306965 
251750258J 



4151708i 
-149835 
-10510 



994043048211675 
-30629467750426^ 
3603111185211j} 
-98042635537 



11784625f 
-9271523*- 
328405 
13S74 
-S42 



48323690206458975 
-159S349963491021i 
151S03540464029f 
-4951488221400j 



23431764751- 
43741046i 

179017654 



Table 3. Coefficients 0,2, s,j and &2,sj in the extrapolation formula (|ll| ) for Z2. 



s/j 1 2 ;i 4 5 6 7 8 9 10 1 1 12 




199008 


35859664 


3105873216 


366981638400 


43962925841280 


5412646450270080 


689859653042488320 


91351020237486105600 


12589478624473041715200 


1806730592586035116492800 


269973775259474769527808000 


1 12 724 


3872 


-3918624 


-664081632 


88392730560 


-11084244867072 


-1384694987443488 


176331199274687232 


-23140865178023786496 


-3146715137132668707840 


-444530245545400234049280 


-65310818282764528381701120 


2 -20 -444 


-34224 


-1984608 


-100703480 


- 3423052408 


152494743744 


58829587491792 


10515670803566400 


1629900454944139200 


243735746832642107520 


36500387580624231020736 


5565229584050566554128640 


3 16 -48 


1788 


325744 


27280764 


2042841660 


148133624112 


10422978961984 


671732365773200 


31014499461745392 


-1105330916217243648 


-683352705387601823616 


-155665991250042317339904 


4 2 62 


2056 


49792 


E47B80 


-62926964 


- 9069431144 


-924884356296 


-87052S467572S0 


-79808557490692SO 


- 71 E3882S 178813 1360 


-6 082 107154338879 2 75 2 


- 4414115889872814761376 


5 2-4 


-286 


-11426 


-409240 


-14714408 


-466639676 


-4246070624 


1552859277500 


259271668545748 


32421044232235536 


3730S95758411050000 


413972960682773319120 


6 -2 




754 


32534 


1427150 


66602412 


3204308732 


144793941316 


4584284710972 


-173506172385408 


-646467S6919141120 


-10662860072358795800 


7 


2 


-14 


-1558 


-67700 


-3405654 


- 188754426 


-11056264376 


-654356019664 


-36883161160012 


- 1658325705468872 


1972409671354252 


S 




-2 


22 


2784 


116304 


6427038 


414871986 


28039294862 


1943327703680 


135094260867632 


8967575039268340 


9 






2 


-32 


-4512 


-174028 


-10163570 


-774000120 


-58982092862 


-4593416253250 


- 365367972622928 


10 










2 


6816 
-58 


232224 
-9764 


13674990 
-277408 


1290331376 
- 15198942 


109505894322 
-1993340778 


9324373410550 
-186906331802 


IS 

13 












-S 


74 

1 


13418 
-92 


290880 
-17834 


11939322 
-248470 


2942846678 
141926 


14 
















-2 


112 


23062 


120410 


15 


















2 


-134 


-29146 


16 
17 






















158 




s/j 1 




3 


4 




6 




8 


9 


iO 


11 


12 





-40 


-1848 


-72579 




- 118290915 


-5121342765 


-235010259540 


-11456914728105 


- 592786074379275 


-32513880609579600 


-1887019647180154987^ 


1 4 


-50 


-1176 


-26112 


-647905 


-16453117^ 


-425308222^ 


-10281313150^ 


-181320757591^ 


2606202716040 


622231248666307| 


565600970776 75104f 


2 -4 -12 


-40 


-344 


-3513 


-44366 


-486457j 


-8785 2 7f 


105956879 


996828993 


-302888778255 


-35962386099941^ 


-2973976984644465}! 


3 4 


14 


-28 


-36 


-3009 


-85659^ 


-1339020 


-70112132f 


-2289397563 


-61537586981| 


-14176450368171 


14046475617195 






-10 


325 


2170 


29241} 


736029 


5619955f 


250213920^ 


8870308461 


182363269816f 


5892713122080 


5 









-945 


-5052} 


-5841} 


-1080203^ 


34632 655^ 


705890175 


1850093234 


713236878450} 


B 
7 
S 

S 








16 


1955 
-38 


4147} 
-3395 
66 
4 


-292766j 
9954 
5278 

"100 


-1058909^ 
131 090 Oi 
-52674 
-;^90 


-217483272} 
5720061| 
-3735338i 
146385 


-5310678979} 
539675633} 
-343039} 
8459448} 


-1036978067} 
13147234615f 
-464689558-5^ 
-61062655} 


ID 














-4 


14Q 


10290 


-321145 


-16469070| 


11 
















4 


J86 


-13310 


615230 


12 


















.4 


238 


16555 
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It is clear from equation (|Tl|) that the — s — 2 terms available from the 
various sequences for the correction terms are not sufficient to determine all the 
2(s + 72 + 2) + 1 unknown coefficients in the extrapolation formulae for large s. However, 
from tables we immediately see that the leading coeffients, an,s,s+n+3 and bn,s,s+n, 
in the extrapolation formulae are alternating in sign but otherwise constant, e.g., 
cin,s,s+n+3 = (— 1)*^"'2 and bn,s,s+n = — (— 1)'*^"'4. In general we have found that the 
leading coeffients are expressible as polynomials in s, 



■1) 



{2s)\ 



s+n 2k+l 



and 



■11 



2k 



''n,s,s+n- 



(12) 



(13) 



j=0 



where / = max(A; — 3, 0). The coefficients of these polynomials are listed in table ^ and 
table |. 

This time we note that the 2 or 3 leading coefficients are independent of n. And 
indeed we find that /3n,fc,2fc-j/3'^ is a polynomial in k of order 2j + 1. In particular we 
have, 

f3n,k,2k = -S\k+l), 

/3n,fc,2fc-i = 3'=(278P + 99k^ - 179A;)/27, 

3'=(77284A;5 - 233636A;^ + 145247P + 233636P - 222531A;)/1458 n = 

3^'(77284A;^ - 233636A;^ + 148487P + 233636^ - 225771A;)/1458 n = 1 

3'=(77284A;5 - 233636A;^ + 151727^^ + 231692P - 230955A;)/1458 n = 2, 



f3n 



fc,2fc-2 



and 



Pn. 



fc,2fc-3 



3^(107424760fc^ - 981604100fc^ + 3689847622A;^ - 5987330165A;^ 
+1103673490A;3 + 6968934265^:2 - 4900945872fc)/590490 

3'=(107424760fc^ - 981604100/^^ + 3703358422A;5 - 6035726045A;^ 
+ 1135273210A;3 + 7017330145A;2 - 4946056392A;)/590490 

3'=(107424760A;^ - 981604100fc^ + 3716869222^^ - 6092228405A;^ 
+1180199050fc3 + 7073832505P - 5004493032A;)/590490 



n 



n 



n 



2. 



So when calculating the extrapolation formulae (|Tl]) we ffist used the sequences for 
the correction terms to predict as many polynomials as possible. When we ran out 
of terms we then predicted as many of the leading coefficients from (|T^ and (|TB|) as 
possible. This in turn allowed us to find more extrapolation formulae, which we could 
use (together with the formulae for Pn,k,2k-j) to find more of the formulas for the leading 



Table 4. The coefficients cxn,k,j in the extrapolation formulae (jl. 



k/j 
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2 


3 


4 


5 












•J 







2 


12 


96 


25920 








1 




2 


-1468 


510840 


1300549152 


6977845103040 


2 




2 


-558 


-888576 - 


■2703028424 


-15961597601040 


3 






-428 


600150 


1803808536 


12990760770024 


4 






6 


38850 


-409368974 


-4669644914724 


5 








-7950 


-18927216 


667412568090 


6 








126 


5412484 


519181770 


7 










-224616 


-5708659068 


8 










6354 


428756328 


9 










-96 


-15678630 


10 












484770 


11 












-10656 


12 












96 










aii.k 


.J 







2 


10 


312 


7200 


967680 





1 




2 


-1192 


566700 


1554817152 


8201223021600 


2 




2 


-522 


-1082526 - 


■3148467656 


-18452430302760 


3 






-428 


592890 


2028534576 


14630011358544 


4 






6 


39480 


-411219494 


-5057308922574 


5 








-7950 


-20088936 


678397444650 


6 








126 


5432476 


1825419990 


7 










-224616 


-5791418388 


8 










6354 


430798878 


9 










-96 


-15704550 


10 












484770 


11 












-10656 


12 












96 










012. k 


.J 







2 


4 


384 


-14400 


483840 





1 




2 


-940 


679560 


1809266592 


9448276380480 


2 




2 


-510 


-1261176 - 


■3582329672 


-20938915839600 


3 






-428 


598470 


2221265016 


16212179478024 


4 






6 


39930 


-415205294 


-5412395679564 


5 








-7950 


-20805456 


691017125850 


6 








126 


5445412 


2648223930 


7 










-224616 


-5853962268 


8 










6354 


432269568 


9 










-96 


-15721830 


10 












484770 


11 












-10656 


12 












96 
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Table 5. The cocfBcients (3n,k,j in the extrapolation formulae (13) 



k/j 


1 


2 


3 


4 


5 


6 


7 


8 


0O.k,j 





-1 























1 


22 


1658 


334656 


153601904 


109900869120 


119835268440320 


179701691558215680 


358810347942675191808 


2 


-6 


-2385 


-649384 


-348301926§ 


-276393269472 


-325054574102944 


-517144503241342976 


-1083429067670151410841 1 


3 




754 


403500 


284418460 


266027954552 


3501885380126875 


606338380910234368 


1357922016436272064896 


4 




-27 


-96524 


-108269551 § 


-129814389840 


-200292896954500 


-388885779063580337 | 


-950703086097155610426 § 


5 






7860 


20243976 


35294435450 


67835150491952 1 


153240568087470602 § 


419124576286737422400 


6 






-108 


-1748436 § 


-5461778850 


-14237073600841 


-39235294734473786 § 


-123701080180932859946 ^ 


7 








55980 


465582548 


1869074318751^ 


6694588679545816 


25258034274069075192 


8 








-405 


-19729020 


-150316130097 


-764582679410893 J 


-3623005587850488121 


9 










326970 


6996149842 1 


57616149115192 


365799796385163936 


10 










-1458 


-168308595 


-2763401680710 § 


-25727820621001866^ 


11 












1688526 


78924402909^ 


1230228773135952 


12 












-5103 


-1196493480 


-38326543230336 § 


13 














8015112 


725790362976 


14 














-17496 


-7483535892 


15 
















35779320 


16 
















-59049 






-1 























1 


22 


1874 


380544 


173844464 


124703585280 


135892573287680 


203891746284011520 


407060480146447669248 


2 


-6 


-2505 


-717448 


-387796806 § 


-309968480352 


-365371396244512 


-582626355682555904 


-1221984137726510908569 1 


3 




754 


427116 


309285900 


293508255992 


388934364085487^ 


676720632830112000 


1519953696088914094464 


4 




-27 


-97964 


-114431071 § 


-140364933360 


-219200100776580 


-429068174740539313 § 


-1054354166560941235770 § 


5 






7860 


20801176 


37316508410 


72997883541272 1 


166848798285151882 § 


459887247973724918464 


6 






-108 


-1759236 § 


-5648008290 


-15047082876273 


-42102339808215418 1 


-134133831039985827818^ 


7 








55980 


472540628 


1940426816751^ 


7075210989955416 


27042955126506979800 


8 








-405 


-19793820 


-153559531937 


-7959363974128773 


-3828497810715795865 


9 










326970 


7058663482 1 


59143290605752 


381544301457679104 


10 










-1458 


-168648795 


-2803243801350 § 


-26506146752612586^ 


11 












1688526 


79382817117^ 


1253639579258096 


12 












-5103 


-1198126440 


-38712986171808§ 


13 














8015112 


728717117184 


14 














-17496 


-7490884212 


15 
















35779320 


16 
















-59049 






-1 8 























1 


22 


1826 


399904 


188173424 


136763942784 


150135943911680 


226399980054487040 


453481872936027199488 


2 


-6 


-2553 


-747352 


-415492486 § 


-337028036672 


-400806757117984 


-643091153929042944 


-1354435396148323267737§ 


3 




754 


438668 


326540300 


315331137752 


422565255743151^ 


741047725232481564| 


1673542368950203659648 


4 




-27 


-98972 


-11871011l| 


-148615548440 


-235371247103300 


-465345164498565169 1 


-1151609375547099787322 1 


5 






7860 


21201176 


38883012026 


77347956643752 1 


178970258611925237^ 


497699046774849639872 


6 






-108 


-1767876§ 


-5793139250 


-15721316702945 


-44622273568721818 § 


-143694733990288250474 ^ 


7 








55980 


478155188 


1999479461647^ 


7405891926422253^ 


28659202377280586328 


8 








-405 


-19848900 


-156259593217 


-822961941093677^ 


-4012551877454570521 


9 










326970 


7112021962| 


60457289130459 i 


395524266443898240 


10 










-1458 


-168948171 


-2837772622950 1 


-27193733020204138 ^ 


11 












1688526 


79787198493^ 


1274324332887664 


12 












-5103 


-1199596104 


-39056845236000 § 


13 














8015112 


731355327360 


14 














-17496 


-7497602676 


15 
















35779320 


16 
















-59049 



9 



coefficients an,s,s+n+3-k and bn,s,s+n-k- We repeat this until the process stopped with 
the extrapolation formulae listed above. 

From Zn,33{u) we extended the series for Zq and Zi to 0{u^^^), while the series for 
Z2 was extended to 0{u^^'^). The resulting new low-temperature series terms are listed 
in table |^. The series terms for n < 80 can be found in The full series is also 
available by electronic mail or via the world wide web (see end of article for details). 





Table 6. New low-temperature series terms for the 


square lattice spin-1 Ising 




magnetisation, M(u), 


susceptibility, xC'"); ^-iid specific heat, Cv{u). 


n 


M(u) 


x(") 


C„(m) 


80 


-54894921926791871723909 


1287288269903730631946751 


14980548594400026006446720 


81 


126297750080690225982572 


-2932618501022211032818300 


-35701518428258787100072584 


82 


-92489409777625802742528 


2052947118516396940212072 


27948358015093163176475128 


83 


-292140800078967381434028 


7265393839008482331992336 


80787294138478683214133578 


84 


1131329573810488998686811 


-27533569403846138701663365 


-328084218297769469375116200 


85 


-1662612103740713574604884 


39920035526199174994626036 


500413572618503616248896480 


86 


-492613608448080934983288 


14931369412552677008749932 


113883610899578636803467152 


87 


8211147564410929129324192 


-209438337886232625955934292 


-2436189367728965310422489418 


88 


-18878709288285563929234997 


476343631011307763277303031 


5796339184388631113650158304 


89 


13785042470417967980505100 


-331962548955170787548255312 


-4521670273968534171069307750 


90 


43860793960742590039383898 


-1182970662366618284903971030 


-13149275656695239669180483520 


91 


-169598161959998919874255236 


4472491833451466220138614096 


53268752093378620751877315482 


92 


249007028837325086293670283 


-6471380268030636441237200097 


-81105273455462683681912082728 


93 


74990743692109664334064752 


-2453067840623244506560874448 


-18811719845845853329570528542 


94 


-1234426063083531162682558560 


34026585140776880977075866240 


395621948197278485674431575608 


95 


2835867624597373150874747480 


-77254798018355427358151075836 


-939772638538221206856049448380 


96 


-2064573852622364936737424098 


53613918339192473420875219089 


730560439047179326970330734464 


97 


-6616142800429142888692342768 


192358286613090067676905433984 


2137403950614601807223987280454 


98 


25541973318840653657077177270 


-725486237294161236962098624258 


-8637803751685796115480314923276 


99 


-37462784919716975353690569292 


1047757642093298413301695350332 


13128860187862258100931770296284 


100 


-11463726370314369427304114523 


402487643299242910269563033131 


3102521140136519939724246308400 


101 


186381597774899493113385553664 


-5520950047108379701377201734452 


-64169713991043524782724303040306 


102 


-427804305737711686514509242390 


12514441143336335641931499176014 


152189962317529977423571879172796 


103 


310497065688338308368382361752 


-8648192214593796849228047528228 


-117900563896358034274507714039430 


104 


1002057377876891572217899261401 


-31239827769954401229313918082438 


-347042122337706823781434388021392 


105 


-3862092565977953294672165562608 


117548646091803611498084800222172 


1399138594951555677686779993658460 


106 


5658355584526190756409048411218 


-169446914107751509953285066975278 


-2122963342745826385809572352025484 


107 


1758697226344654211271247297588 


-65942014127753928407511143000684 


-510975035983101017637996662363534 


108 


-28247139203735773793168037850532 


894839314320840238617134121410478 


10397649426241825241936709101268192 


109 


64775215576079739888972057107032 


-2025063464948777359185732816409016 


-24621560654630485557870861813977290 


110 


-46866364052133130252606723536110 


1393635387942882646442005194185966 


19008526767153717074391591041349580 


111 


-152306338289478778693447758360664 


5068325152947612595253210347568392 


56293290314061858138013078241465190 


112 


586011194798187595036300267855234 


-19027185558928704352441507167871531 


-226417646512838638150966074051275808 


113 


-857572414389308013646509973160724 




342972424908794790373332519783688916 



4. Analysis of the series 

We analysed the series using the same methods as in our previous paper [Q] to which 
we refer the reader for details. Here we will give only a short summary of the results 
including improved estimates for the critical point and amplitudes. 
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The estimates Uc = 0.5540663(5) for the physical singularity and /5 = 0.12507(2) for 
the critical exponent of the spontaneous magnetisation were obtained from homogeneous 
differential approximants (which are equivalent to Dlog Pade approximants) by 
averaging over [A^, M] approximants with \N — M\ < 1 using at least 100 series terms. 
The figure in parenthesis represents the spread among the approximants (basically one 
standard deviation) and should not be viewed as a measure of the true error as they 
cannot include possible systematic sources of error. From these estimates it is clear 
that f3 = 1/8 as expected. However, the estimates converge very slowly towards this 
value and even with a series as long as the present 114 terms the estimates have not yet 
settled down to their true value and there is a slight downwards drift in the estimates 
for both Uc and (3. Analysis of the susceptibility and specific heat series yield exponent 
estimates fully in agreement with the expectations that 7' = 7/4 and a' = 0. By using 
our knowledge of the exact values of the critical exponents and assuming that close to 
Uc the estimates for the exponents depend linearly on the estimates for the critical point 
(inspection of the various approximants clearly supports this assumption) we are led to 
the improved estimate for the critical point Uc = 0.5540653(5). 

We calculated the critical amplitudes using two different methods, both of which 
are very simple and easy to implement. In the first method, we note that if f{u) ~ 
A(l — u/uc)^^, then it follows that (uc — u)f^^'^\u=uc ~ A^^'^Uc- So we simply form the 
series for g{u) = {uc — u)f^^^ and evaluate Pade approximants to this series at Uc- The 
result is just A^/^Uc- This procedure works well for the magnetisation and susceptibility 
series (it obviously cannot be used to analyse the specific heat series) and yields the 
estimates Am = 1.20840(5) and A^ = 0.06172(4) where the error bar primarily reflects 
the uncertainty due to the estimate of Uc- For the specific heat series two different 
approaches have been used. In the first approach we look at the derivative of the 
specific heat series for which the above method should work with A = 1. This yields the 
estimate Ac = 22.3(1). In the second approach we start from f{u) ~ Aln(l — u/uc) 
and form the series g{u) = exp(— /(«)) which has a singularity at Uc with exponent 
A. One virtue of this approach is that no prior estimate of Uc is needed. However, 
the spread among estimates from different approximants is quite substantial, though 
the amplitude estimate is consistent with that listed above. Biasing the estimates at 
Uc yields Ac = 22.3(3). In the second method, proposed by Liu and Fisher [^], one 
starts from f{u) ~ A{u){l — u/uc)~'^ + B{u) and then forms the auxiliary function 
g{u) = (1 — u / UcY f {u) ~ A{u) + B{u){l — u/uc)^- Thus the required amplitude is now 
the background term in g{u), which can be obtained from inhomogeneous differential 
approximants [||. This method can also be used to study the specific heat series. One 
now starts from f{u) ~ A{u) ln{l — u/uc) + B{u) and then looks at the auxiliary function 
g{u) = f{u)/ ln(l — ujUf). As before the amplitude can be obtained as the background 
term in g[u). This analysis yields the amplitude estimates listed in table |^. With 
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the exception of the specific heat amphtude there is excellent agreement between these 
estimates and those obtained from the first method. 

Table 7. Estimates for the physical critical amplitudes Am, A^ and Ac from 
the method of Liu and Fisher obtained from inhomogeneous first-order differential 
approximants. L is the degree of the inhomogeneous polynomial. 



L 


Am 




Ac 


5 


1.208284(74) 


0.061712(47) 


20.20(26) 


6 


1.208270(15) 


0.06173(12) 


20.13(22) 


7 


1.208256(80) 


0.06170(20) 


20.23(23) 


8 


1.208266(16) 


0.06164(21) 


20.29(20) 


9 


1.20824(21) 


0.06170(12) 


20.28(32) 


10 


1.208269(35) 


0.061866(60) 


20.31(26) 


15 


1.208272(74) 


0.06165(15) 


20.47(32) 


20 


1.20824(12) 


0.06167(15) 


20.464(75) 


25 


1.208250(47) 


0.06172(11) 


20.448(46) 


30 


1.208236(51) 


0.06174(20) 


20.42(12) 


35 


1.208228(60) 


0.061792(82) 


20.40(20) 


40 


1.208263(70) 


0.06159(20) 


20.45(12) 



Regarding the value of the confluent exponent Ai we have little to add to our 
previous results. Even with a series as long as 114 terms we could not obtain accurate 
estimates for Ai. Again the Baker- Hunter ||^ transformed series of the magnetisation 
favours a value around 1.05 while estimates from the suceptibility series again fall in 
two groups around 1.15 and 1.4, respectively. Using the transformation of Adler et al 
i 

G{u) = \F{u) + (ue - u)dF{u)/du, 

where F{u) is the original series and A the leading critical exponent, yields estimates 
consistent with Ai = 1 for both the magnetisation and susceptibility. 

We find a non-physical singularity closer to the origin than at u±. = 
-0.3019395(5) ± 0.3787735(5) with exponents p = -0.1690(2), 7' = 1.1692(2) and 
a' = 1.1693(3), and a singularity on the negative w-axis at U- = —0.598550(5) with 
exponents equal to those at the physical critical point. Note that our estimate for (3 
at u± has changed substantially from that given in our previous paper. This is mainly 
because we have put greater emphasis on estimates obtained from inhomogeneous first 
and second-order differential approximants. We note that we now have firm evidence 
to show that the scaling law a' + 2/3 + 7' = 2 holds at both the physical as well as the 
non-physical singularities. 
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E-mail or www retrieval of series 



The series for the spin-1 Ising model can be obtained via e-mail by sending a request to 
iwan@maths.mu.oz.au or via the world wide web on [http: / /www.maths. mu.oz.au/~iwan/ 
by following the instructions. 
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